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SUMMARY 


The virell-known methods of thin-airfoll theory have 
been extended to obllqxie or swept-back airfoils of finite 
aspect ratio moving at supersonic speeds. The cases con- 
sidered thus far are symmetrical airfoils at zero lift 
having plan forms bounded by straight lines. Because of 
the conical form of the elementary flow fields the results 
are comparable in simplicity to the results of the two- 
dimensional thln-airfoil theory for subsonic speeds. 

In the case of untapered airfoils swept back behind 
the Mach cone the pressure distribution at the center 
section is similar to that given by the Ackeret theory 
for a straight airfoil. With increasing distance from 
the center section the distribution approaches the form 
given by the subsonic-flow theory. The pressure drag is 
concentrated chiefly at the center section and for long 
v/ings a slight negative drag may appear on outboard sec- 
tions . 


INTRODUCTION 


In reference 1 it was pointed out that the wave drag 
of an infinite cylindrical airfoil disappears when the 
airfoil is yawed to an angle greater than the Mach angle. 
This observation led to the conclusion that the drag of 
a finite airfoil could be greatly reduced by the use of 
sufficient sweepback. With such a swept-back wing the 
wave drag would be associated with departures from the 
ideal two-dimensional flow at the root or tip sections 
and would thus be a -function of the aspect ratio. The 
present report extends the theory of reference 1 to take 
account of these effects. 
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The treatment is based on the theory of small dis- 
turbances in a frictionless compressible fluid. The 
idealized fluid and its equations of motion are identical 
v/ith those employed in acoustics in the theory of sound 
v/avss of small amplitude . The applice.tion of the theory 
is thus limited to bodies having thin cross sections so 
that the velocity of motion imparted to the fluid is small 
relative to the velocity of sound and so that the pres- 
sure disturbances produced are small relative to the 
ambient pressure. 

Ihe adaptation of the sound-vjavo theory to the aero- 
dytiarvics of moving bodies v;as suggested many years ago 
by Prandtl. The theory was apnlied by Ackeret (refer- 
ence 2) to thin airfoils moving at supersonic speed, 
Ackeret ’s treatment is limited, however, to infinitely 
long cylindrical airfoils moving transversely. The 
present theory may bo considered an extension of Ackeret ’s 
theory to take into account wings of finite span and wings 
havirig tapered or swept-back plan forms . In the case of 
swept-back plan forms* the results are markedly different 
from those obtained by tho Ackeret theory and approach 
tho V .lues indicated in references 1 and 5* 

In reference Busemann describes a method for cal- 
culating the supersonic flov/ over bodies v/hich produce a 
conical pressure field. Busemann shows that the flow 
around cones of circular cross sections as well as the 
flov7 around the tip of a rectangular lifting sur-face 
satisfies this condition. Tho fact that a great variet3’’ 
of three-dimensional flows can be constructed by the 
superposition of conical and cn^’llndrical flow fields 
leads to an essential simplification of the airfoil 
theory at supersonic speeds. 

"fhe present treatment differs from Busemann 's in 
that it is further limited to flat bodies, that is, 
bodies vdiich are thin in both longitudinal and transverse 
sections . This additional restriction leads to a much 
simnler mathematical treatment and one which is apu lie able 
to a v/ide variety of airfoil shapes. In a paper presented 

before the Institute of Aeronautical Sciences^, .Allen iU. 
Puckett also treats symmetrical non-lifting bodies. 
Puckett's method makes use of integral expressions corre- 
sponding to tho velocity uotential of pls.nc-source distri- 

bu t i on . _________ 

"llpth Annual Meeting, Mew York, January 29, 19l|.6 . 


2 


IIACA TIT ITo. 1107 


y 

SY 1 -I 30 I/S 

flight velocity 

T'T 

A.* 

Mach number 

S 

coordinates 
point on X-axis 


limit of integration 

dis turbance -ve loci by potential 

U, Vj V/ 

dis turbax.ce -ve loci ty components 

n 

value of u at xj 

u 

value of u for conjugate arrangement 

p 

local pressure 

q 

d^niamic pressvu’e 

p 

density of air 

P Q 

Legendre functions 
source -strength factor 

D 

differential operator 


drag coefficient 

t 

m 

thickness of wing 

slope of lino source (absolute value) 

c 

chord of v/ing 

TIiS OBLIQUi LINS SOURCE 

The 

velov'^ity 

assurrptions of small disturbances and a constant 
of sound throughout the fluid load to the 
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well-known linearized equation for the velocity poten- 
tial </j (see reference 5) 


(l - + 9ljy * 9^zz = ^ 


( 1 ) 


The analysis Is simplified by introducing the coordinates 


Xt = X 


y-j^ = \jlA^ - 1 y 


2 


Zj - - 1 z 


( 2 ) 


Dropping the subscripts from the transformed coordinates 
gi ves 


^xx “ ^yy ” ^zz ^ 


(3) 


According to the thln-airfoll theory the pressures 
on the transformed airfoil are given by 


Ac = 2- = - ~ 


•V V 


Ox 


(z — )-0) 


(^) 


dz 

and the slope of the airfoil surface — is equal to the 
slope of the streamlines near the chord nlane ; that is, 

dz _ w 
dx “ V 



(z— >0) 


(3) 
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The elementary solution of equation (5) foi’ ^ point 
source is 


1 



This solution is directly related to the subsonic poten- 
tial 




1 

+ y^ + 


In the subsonic case the equlpotential surfaces are, hov/ever, 
ellipsoids, whereas in the supersonic case the equipo- 
tential surfaces are hyperboloids limited by the Mach 
cone. (See reference 5 derivation of these ele- 

mentary solutions.) 

Because of the linearity of equation (1) a solution 
may be used to denote one of the velocity components 
rather than the velocity potential. The specification 
of one component in this mariner actually describes the 
whole flow field since the other components may be obtained 
by Integrating the given component to obtain the velocity 
potential and then differentiating the results along the 
desired directions to obtain the desired components. 

This procedure is especially useful in the thln-airfoil 
theory, where the complete velocity field may not be 
required . 

Adopting the foregoing procedure, one may write 


yx2 _ yii _ 

Since u is proportional to the pressure, such a solution 
corresponds to a point source in the pressure field. The 
solution for an oblique line source may be obtained by 
integrating for the effect of a row of point sources along 
the line y = mx. it will be shovm that such a line source 
satisfies the boundary condition for a thin wedge-shape 


5 


.v! '- V/* .‘.X. 
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bodjT'. This solution, as well as other expressions relating 
to oblique airfoils, can be most conveniently expressed 
by referring to the oblique coordinates 


X ' = X - my 


y' = y - mx 
z ' = \ll - in^ z 


(See fig. 1.) It may be shown that if any function 
f(x, y, z) is a solution of 


f 


XX 



0 


then f(x', y' , 2 ') is also a solution. In particular, 
the point-source solution becoraes 



Hence the Integration for the effect of an inclined line 
of sources may be nerformed directly along the oblique 
x’-axls; thus, for ra< 1.0 


u = 


I 

to 


1 d§' 

x’ - 5 ' ) ^ - y ' ^ - 2 ' 2 


= I 


cosh 


X ’ 


^y + z ’ ' 


( 6 ) 
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'.vliere is the position of the last soiATce whose ilach 

cone includes the point (x', y*, a*) and is given by 


^ 1 


= x» 


-f 


.2 


a'2 


Vilisn m an-xroaches 1.0 the source line apnr caches 
coincidence v'ith the Mach cone, corresponding to a trans- 
verse velocity component equal to the velocity of sound. 

For values of m greater then 1.0 the integration 
yields 


t 


u = -licos 


t - + 2 ? 2 


(7) 


It will be seen that in this case I is imaginary 


T}ie vertical velocity near 2=0, v/hich determines 
the shape of the boundary, may be determined by inte- 
Curating u with respect to x and then dif f erenuiating 
the resulting velocity potential v^ith respect to z; 
thus (see appendix) , 


w 


60 " _ 6 
6z 6z 


u dx 

J 


- iTT 


m 


\(t - m 2 


( 8 ) 


if z — >0 and y’ < 0. If y' > 0, w = 0. There is 
thus a discontiriuit:/ in the vertical velocity of the 
streamlines when the^/ cross the line source at = 0. 

For small values of l/m this discontinuity in vertical 
velocity’' agrees with the boundary condition for a simple 
wedge shape having a small wedge angle. (See fig. 2.) 
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Tf the source strength I is held constant and in 
is allowed to aonroach sero, the wedge angle ultimately 
becomes large. At m = 0 the line source actually 
satisfies the boundary condition for the circular cone 
(reference 6 ) , but it is found that the slope of the 
conical boundary does not agree with the slope of the 
streamlines near z = 0 and hence the theory no longer 

in 

holds. The condition j — ?0 thus represents the transi- 
tion from an oblique airfoil to a body of revolution and 
will be avoided in the present analysis by restricting 
the formulas to flat bodies, that is, airfoils that are 
thin in both longitudinal and transverse section. 


AIRFOIL OF VvEDGS SECTION 

Over the wedge section near the plane z = 0, the 
forraula ( 6 ) becomes slmpl^y 


u = I cosh' 


JLL 

!y 


( 9 ) 


where |y' j denotes' the absolute magnitude of y' = y - mx. 
Tiie pressure is thus constant along the radial lines 


— = Constant (10) 

Y ’ 

ana is conveniently represented by the variation along a 
line parallel to the x-axis. Figure £ shows the oblique 
wed^je-shape figure corresponding to a line source with 
m< 1.0. In this case the pressure field is confined to 
the interior of the Mach cone - y£-s2 = x'£-y' 2 _ 2 ;' 2 =Q 

and the theory, unlike the Ackeret theory, indicates a 
stagnation point along the leading edge. (x^ctually, of 
course, the thin-airfoil theory shov/s an infinite velocity 
at such ooints, but this is to be interpreted as a velocity 
of the order of magnitude of the flight velocity V. The 
pressure to be expected along the leading edge is the 
stagnation pressure corresponding to the transverse 
velocity component.) 
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Given 


dz w 


dx 


— , the ¥/edge angle measured in do¥/n- 

V 

stream sections, -the souj^ce strength must vary with m 
according to 


_ V m. dz 
“ ” TT I z Cix 


\r ; 

yj. - m- 


( 11 ) 


(from equation 7)* Daen 


Ap _ 2 dz 
q “ IT dx 


in 


cosh 


-1 


- m2 1'^’ I 


U2) 


If m exceeds 1.0, the leading edge of the airfoil 
vrlll lie outside the Mach cone. In this case 


A _2 _ 2 dz 

q 


m 


ym‘- - 1 


cos 


-1 


X' 


I I O O 

+ Z»2 


( 15 ) 


In the region bet\:een the leading edge and the Mach cone 
cos~- ~ — is constant and equal to rr; hence the pres- 


!y’ I 

sure in this region is constant, that is. 


Ap 

q 


= 2 


as m 


dx 


\1 ^ 

yrn- •• 


(lli.) 


Figure 5 illustrates this result. 
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If m — >co a semi -infinite airfoil with its leading? 
edge at right angles to the dii’ection of flight is 
obtained; here 


•X - my -y 


^(y - mx)^ + (1 - m^)z^ Vx^ - 


and ~ = 2— wherever v > yx^ - z^. This value 
q dx - If 

with the Ackeret theory. 


(15) 


agrees 


AIRFOILS BOUNDED BY PLANE SURFACES 


The distribution of pressure over symmetrical air- 
foils bounded by plane surfaces can be obtained by super- 
imposing the pressure fields for several line sources and 
sinks. This superposition is greatly simplified by the 
conical form of the pressure field for each single line 
source. Because of this form the whole distribution in 
the plane z = 0 is, in effect, represented by a single 
curve. If the velocity field for a line source beginning 
at the origin (equation (6)) is denoted by u and that 

beginning at x = -1 is denoted by u_j_, and so forth, 
the sum 


u-i - u+i 

represents the velocity over a plate of uniform thickness 
having a beveled leading edge of constant width. (See 
fig. q.) Similarly 


u_x - 2u + u+ 3 _ 

represents the pressure field for an airfoil having 
diamond-shape cross sections. 

The superposition required for several sources or 
sinVcs can be accomplished by manipulation of a single 
curve if it is remembered that u is a function of the 
ratio x/y. Figure ij. illustrates this process for a 
source and a sink. In terms of the ratio x/y the 
separation of source and sink and hence the scale of the 
chord length continually diminishes with increasing 
distance from the root section. 


10 
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At Ici-rge distances from the vertex (x'— jcx>) the 

e-Toi'ess ion (for in < 1.0) 

,_i x'+l _-| x’-l 

u_j_ - cc cosh ~ - cosh (16) 

jy'-m) jy'+mj 

is found, to approach the value 



where Q,q is the Legendre function (see reference 6). 

In the thin-airfoil theory for subsonic speeds it 
can be shown ths.t if 

w ?n(x) 

(12) 
( 19 ) 

since Neumann's formula (reference 7> P* II 6 ) 

/H-1 

1 Pn(&) 

Qn = p (20) 

may be interpreted as the integration for the velocity 
distribution due to an array of sources of strength 

wdg = Pn(§) dS 


then 


OC \T 


oz 

6x 


u GO i4^(x) 
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along the chond of airfoil. The expression Q.^{j'An) 
of equation (I7) tiim represents the subsonic pressure 
distribution over tlie beveled edge. ^ 

At the root section iy ~ 0 ) only the forward source 
need be considered, since the airfoil STitface is ahead oi 
the Mach cone originating at the rear so\u?ce . Here 


^ GOsh 


.X X -i- 1 - my 
y - m (x + 1) 


G& cesh 




( 21 ) 


and the pressure ove’r th© root section is thus constant, 
as given by the Aekeret theory', bvit is altered in magnitude 
by the obliquity. 

The oblique wing lying behind the Mach lines thus 
shows the Acker et type of pressure distribution over the 
foremost section and a progressive change along the span 
from this distribution to the subsonic type of distri- 
bution. Since the subsonic type of distribution shows 
no pressure drag, there is a continuous falling off of 
the pressure drag with increasing distance from the root 
section. The pressure drag of the oblique wing thus arises 
chiofly on the foremost section, and it fo 3 . 1 ows that the 
drag coefficient of the wing as a whole diminishes with 
increasing aspect ratio. It vrill be sh ovm sul^sequently 

^Similarly, if is taken as the chordwise 

distribution of vortlcity. 


uon ?„(x) 


w cc Q^(x) 


The first of this series of airfoils is. the cambor shape 
cuT’ved to support a \anifoi*m load. 
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that the effect of cutting the wing off along a line 
j = Constant to oroduce a dovmstream tip causes a reduc- 
tion of the pressure drag on the adjacent sections, and 
if the aspect ratio is sufficiently high, the pressure 
drag in the region of the downstream tip may actually be 
negative , 

If the wing lies ahead of the Mach lines (ra > 1.0) 
the Ackeret type of pressure distribution occurs and a 
pressure drag arises over the whole length. In this case 
both u and w are constant over the beveled part at a 
distance from the origin. 

The treatment thus far applies to serai -Infinite 
cylindrical v/ings having root sections near the origin. 

A complete swept-back wing may be obtained by the addi- 
tion of a symmetrical or conjugate arrangement of source 
lines below the x-axis. Values_of u for this conjugate 
arrangement may be denoted by u. Figures 2 and 5 show u 
for a single inclined source and figure 5 shows calculated 
pressure distributions at several sections along the span 
for a com.plete swept-back airfoil having beveled sections. 
The addition of the conjugate source lines doubles the 
pressure at the root section, but this Interference effect 
falls off rapidly along the span. It is noted that, as 
in figure I 4 ., the most significant change in pressure 
distribution occurs along the expansion wave originating 
at the trailing edge of the root section. Figure 6 shows 
the variation in pressure drag along the span for this 
airfoil obtained by integrating the chordwlse components 
of pressure at the different sections. 

The addition of a reversed source-sink distribution 
having Its origin displaced to a point O 2 (see fig. 7) 

will show the effect of cutting the wing off in a direc- 
tion parallel to the direction of flight, it v/ill be 
evident that the effect of such a 1:lp is characterized 
by the subtraction of the curves u and is limited to 
the area lying within the Mach cone which originates at 
the tip. It is interesting to note that pressure distri- 
butions of the Ackeret type, except reversed in sign, 
are added near the tip; hence, cutting the tip off in 
this manner reduces the drag of adjacent sections. 

Figure 8 shows the pressure distributions over a 
rectangular airfoil having a leading edge at right angles 
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to the flow. In the triangular area ahead of the Mach 
cones originating at the tips the pressure is constant, 
as ^Iven by the Ackeret theory, whereas behind these Mach 
cones the pressure drops sharply. 


AIRFOIL OP BICONVEX SECTIONS 


Curved surfaces require a continuous distribution 
of sources and sinks alined with the generators of the 
surface. Each elementary source-line causes an infini- 
tesimal change in direction of the surface and hence the 
slope at any point may be obtained by adding up the 
effects of all sources ahead of that point. Thus 


or 


dz 

dx 


nx 

TT 

V 

iJxo 




( 22 ) 


d^z _ TT \h. - m^ dl 
dx*^ V m ox 


For airfoils of constant chord, m will be a constant 
and the integrations can be performed without difficulty. 
The simplest case is that of constant curvature, which 
leads to profiles formed from clrciTlar arcs. 


In order to obtain a biconvex profile it is necessary 
to Introduce finite sources of strength sufficient to for-m 
the desired angle of intersection of the arcs at the 
leading and trailing edges, together with a uniform dis- 
tribution of sinks along the chord line between the two 
sources. These profiles thus require a uniform distri- 
bution of sources or sinks, which may be obtained by 
Integrating the elementary solution for the line source 


( equation 


by 



( 6 )). 


The resu.lt Ing solution ma^/ be denoted 


lit 


and is, for m < 1, 
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nx 


"U= I 
D 


Cosh 


-1 S - my 


mci 




(■ 


\fe_ 


m*- 


m 


•y 


cosh' 


-1 


l^i 


m 


cosh 



(23) 




jyl 


form 


Inasmuch as the elementary solution u 

/ X \ 

f/ — ] the Integrated solution appears 


is of the 
in the form 


1 

— u = y X 

D 



and will be conveniently represented by a curve typical 
of all snanvvlse stations, namely. 



For a closed profile Intersecting the X-axis at the 
points ±1 there Is obtained 



^'-1 + ^+1 




( 21 ,) 


This superposition may be 
transposing and adding the 
as shovm in figure 9* 


accomplished conveniently by 
typical curves u and — u, 

yi> 


It will be found that if m. is less than 1.0 the 
velocity distribution approaches, with increasing distance 
from the root section, the form given by the subsonic-flow 
theory for an airfoil of biconvex section, that is. 


w 



u cr. 



( 25 ) 
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A 1 r. ths root soction, hov/evSi*, th© rorn is siTsiply theft 
given by the Ackeret theory for a straight airfop, 
although the values are reduced in raagnitude oy the 

m _i 1 

factor -yemm: cosh — • 

\Jl - ni~ 

The pressure distribution and the variation of drag 
along the^ span for the syr.imetrical biconvex wing are snevm 
in figVLres 10 and 11 . 


CONICAL STJRPAGSS 


For tapered airfoils both m and I v/ill be 
functions of (See equation ( 2 '}-).) It is easily seen 

that closed siu'facss can be obtained only on the condi- 
tion that the line sources have a coranon point of inter- 
section, as in figure 7. If this point is denoted by 

9 if o 


m = 





The surface obtained is 
through the fixed point 
surface . 


one generated by a line passing 
'-oj^No hence is a conical 


The pressure ever the tapered airfoil requires the 
integration of 


u = 


-1 - my ^ 

jy - m (x - &) 1 dg 


cosh 


T v 


dl 




v±iere 


Is the location of the vertex of the airfoil and 


dl _ V m dA 

as dg2 


In conclusion it shoiild be noted that the press^ores 
have been derived for an airfoil transformed according 
to equations (2). The pressures at corresponding points 
of the original airfoil are to be obtained by dividing 

by H^- - 1. 


Langley Memoria.1 Aeronautical Laboratory 
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APPENDIX 


EVALUATION OP INTEGRAL OP SX^ATION (8) 


For rc < 1.0 the disturbance is zero outside the 
Mach cone and the range of integration should be extended 


;nlY from x 


'I = ^ 


- + Z-- to X, that is. 


u dx = 

t;-<» \j 


IX 


cosh 


hV. “ 1 


\^2 +22 \ fy »2 + 2 , t 2 


dx (Al) 


(for unit soixrce strength). F’..irtherniore, 



u dx = 




( 


du 

bz 


dx 


7-2 + z2 


(A2) 


since the Integrand Is zero at the lo\';er limit. 


iOW 


<=• ^-1 
— cosh 


'T t 


-X • s * \/l - m2 




y'^ + z'2 (y,2 + J., 2 ) - y 


(A5) 


I 2 _ s t c; 


and hence the Integral 


OX. 


w 


-X ' z ’ Vl - 


d; 


\/y2 + 2 2 + z'-) Vx>2 . y.2 . 2*2 


(Ak) 


must be evaluated. 
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First it is noted that the integral vanishes 
with z except in the neighborhood of the Mach cone 



line source 
so that 


- z'2 = 0 ) 

(y ' = 0) . 


and in the neighborhood of tl;e 
Near the Mach cone y’- + z’2-^ 


t2 

9 



-•X ' z ' dx 



X’ 


-z dx 



(A5) 


Since the latter Integral approaches zero v/ith z, there 
is no contribution to equation (Aiq) in the region of the 
Mach cone. On the other hand, near the line source y' — ^0 

and ^x ' ^ - y ' ^ - z'^ >x ' ; hence, as z’ — >0, 





+ z 


-X ' z ’ dx 

2) 


o 




y 

— + Constant 
z ^ 


(A6) 


The value of the integral changes from 0 to tr in crossing 
over the line source at y’ = 0 and is positive or negative 
depending on whethei- z' approaches zero from the posi- 
tive or negative side of the xy plane. Hence 


w = ±— 'll - m~ 
m 


If m is greater than 1.0 


u 


cos 


-1 


X' 



(A?) 


19 


NACA TN No 


1107 


and the flov/ disturbance extends outside the Mach cone 
to a region bounded by plane waves extending from the line 
source and tangent to the Mach cone. (See fig. 12.) llae 
equation of these planes can be easily shown to be 

jtc: ^ - Q. hence for ra > 1.0 the lovrcr Ihnit of 

intx'gration is given by 

yi2 2»2 = 0 


or 


^ 7 , 


X. 





Then 


(A8) 


aX 




A 

6 7. 


u do: 




A I 

oz } 

.i 


iOS 


1 


X 


I 



» 2 


dx 


(A9) 


In this case u does not go to zero at the lov\rer limit 
but is equal to tt. In all other regions, hov/over, 
the integral approaches zero uniformly with z as in the 
preceding case ; hen c 0 


w = 


AX 


6z 


u dx = u 


0x7, 


I 


- + 


TT 

m 




- 1 


(AlC) 


as before. 
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Fig. 1 



Figure I.- Oblique coordinotes. 
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Fig. 2 
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Fig. 3 



Figure 3. — Pressure field for oblique 
wedge where m>io ^ 
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Fif. 4 
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Figure 5” Variation of pressure distribution along span 
of swept— back wing. m*tan 30®. 
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Figure 6.— Variotion of drag coefficient with distance from 
root section for swept~back wing. Wedge section. 

M = l.4. 
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Figure 7. — Addition of reversed source-sink 
distribution to produce tip. 
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Figure 10.- Pressure distribution at 
different points along span. Biconvex wing. * 
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Distance from root section in !/2 chord lengths 

Figure II.— Variation of drag coefficient with distance 
along span for wings of biconvex section. M = 1 . 4 , 

4- =10 percent . * national advisory 
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